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Modelling the Broadband Spectra of MCG-6-30-15 with a 
Relativistic Reflection Model 
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ABSTRACT 

The spectrum of the bright Seyfert I galaxy MCG-6-30-15 shows a broad Fe Ka emission 
line, generally interpreted to originate very close to the central black hole which must then 
have a high spin. The observed X-ray variability is driven by a powerlaw component, with 
little variability found in the iron line energy band. The disconnection of continuum and line 
variability may be a consequence of strong gravitational light bending. The X-ray spectrum 
of MCG-6-30-15 is however complex, strongly modified at soft X-ray energies by warm ab- 
sorbers, which some workers have extended to build an absorption-dominated model of the 
whole source behaviour The absorption model interprets the whole X-ray spectrum without 
a relativistic reflection component and attributes most variability to the warm absorbers. We 
re-examine the XMM-Newton, Chandra and Suzaku data taken in 2001, 2004 and 2006, re- 
spectively, and construct a model consisting of several warm absorbers confirmed by spectral 
analysis, together with a relativistically blurred reflection component that explains both the 
soft and hard excesses as well. The model works well on data from all epochs, demonstrating 
that the reflection model provides a consistent interpretation of the broadband spectrum of 
MCG-6-30-15. 
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1 INTRODUCTION 

The X-ray reflection spectrum of Active Galactic Nuclei (AGN), 
included fluorescent lines, the Compton hump and other emission 
in the low-energy band, is expected from a cool accretion disc illu- 
minated by a power-law continuum iRoss&Fabianl 19931) . As cold 
gas absorbs high energy photons, inner-shell electrons are knocked 
out and, followed by atomic transitions, lead to emission lines. The 
iron Ka line which lies at or above 6.4 keV is usually the strongest. 
This fluorescent line is intrinsically narrow, but can appear as a 
broad line because of Doppler effect and gravitational redshift ef- 
fects dFabian et al.ll989ll200(]|.lRevnolds & Nowak2003.) . A broad 
iron line was first discovered by the Advanced Satellite for Cosmol- 
ogy and Astrophysics (ASCA) in the X-ray s pectrum of the Seyfert 
I galaxy MCG-6-30-15 (Tana ka et alJl 19951) . Since then MCG-6- 
30-15 has become one of the most studied AGNs. The broad and 
skewed iron line was o bserved in later obser vations of BeppoSAX 
I Guainazzi et alj|l999l) . Chandr a jLee et al.ll20 02'). XMM-Newton 
I Fabian et al.ll2002l) and Suzaku jMiniutfi et ^200711 . 

The soft X-ray spectrum of MCG-6-30-15 has complex struc - 
tures usually attributed to warm a bsorbers (e.g. lOtani et 311119961) . 
iBranduardi-Ravmont et al.l ( 1200 ll) tried however to explain the soft 
features by the same mechanism used to interpret the broad iron 
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line. They described the 0.7 keV drop as the blue wing of a broad O 
vm emission line, and the smaller ~0.39 keV and ~0.54 keV drops 
as C VI and N vil emissio n lines, respectively. Re flection modelling 
of the soft X-ray lines (B allantvne et al.|[2002h failed to confirm 
the large equivalent widths required by the data. iLee et al] ( l200lh 
presented evidence of dusty warm absorbers in the Chandra High 
Energy Transmission Grating Spectrometer (HETGS) observation. 
The 0.7 keV drop can be interpreted as a combination of resonance 
absorption lines redward of the O vil edge and an Fe I L-shell ab- 
sorption edge, caused by dust grain iron oxides or silicates. The 
presence of du s t is in agreement wi t h the optical observations of 
[Reynolds et all h997h . 'Turner e t al.l ( l2003l) examined the Reflec- 
tion Grating Spectrometer (RGS) spectra from XMM-Newton ob- 
servation and confirmed the presence of dusty warm absorbers . The 
total X-ray spectrum of MCG-6-30-15 could hence be illustrated 
by a power-law continuum, the reflected component from the cool 
disc, and absorptions by warm absorbers. 

An initially puzzling lack of line variations during rapid 
continuum variat ions (Fabian et al. 2002; Matsumoto et al. 2003|) 
aroused interest. iMiniutti & FabianI 120041) showed that if the illu- 
minating source moves within the innermost area around the black 
hole where there are strong general relativistic light-bending ef- 
fects, then the flux of the continuum varies but the iron line could 
stay unchanged. This gravitational light-bending model inteiprets 
the lack of variability of the iron line, despite the observed variabil- 
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ity of the continuum, as due to changes in the location of the power- 
law component. Further studies of this model have been report ed by 
Suebsuwong et al. (2006), Nie dzwiecki & Mivakawal ^20101) , and 
Niedzwiecki & Zvcki t2008i) . The timing analysis of the Suzaku 
observation in 2006 lat er showed consistenc y with the gravitational 
light-bending model jMiniutti et all |2007|) . Ilnoue & Matsumotol 
(2003), on the other hand, suggested a different view with an 
absorption-dominated model which could explain the X-ray spec- 
tra of MCG-6-30- 15 as well. The warm absorbers might imitate the 
"red wing" structure seen in the iron line energy band, and varia- 
tions of continuum are hence due to warm absorbers but not of the 
power-law component. Miller et al., (.2008, ,2009.) further extended 
the absorption model to build a complex "3-1-2" model which con- 
sists of five absorbing zones, two of which are 'partial-covering'. 
With 3 percent systematic errors included, they claim that this 
model reproduces the variability and XMM spectra of MCG-6-30- 
15. The partial covering component responsible for the red wing 
shape produces no other observed featur es at CCD resolution. Us- 
ing Suzaku data, iMivakawa et alj ( |2009| ) also favour this view but 
reduce the need for five warm absorbers to two. 

Apart from the variability disconnection, it has been claimed 
that the strong emission found at high energies (above lOkeV), 
known as the "hard excess", cannot be explained with a simple disc 
reflection scenario. Miller et all (^2009') concluded that high-column 
partial-covering absorption dominates the hard excess in MCG-6- 
30-15, as seen in N GC 1365 (Risaliti et al. 2009), IH 0419-577 
dTumer et alll2009t) and PDS 456 tReeves et al.ii2009,) . The spec- 
tra of the last three objects have however since been shown to be 
well described by reflection ( Walto n et al.ll2O10h . although Comp- 
ton thick, partial-covering absorbers could not be completely ruled 
out. 

We re-examine the full X-ray data of MCG-6-30- 15 taken 
after 2001 by XMM-Newton, Suzaku and Chandra. A reflection 
model with absorbing zones was applied to these observations to 
test the robustness of the model. We start by investigating the differ- 
ence spectra of the three datasets and then study the full real spec- 
tra. The grating data of Chandra, the energy resolution of which 
is better than other instruments, were used to study the effects 
of warm absorbers and determine their shape. A model with four 
absorbers works well on Chandra and Suzaku spectra, while the 
spectrum of XMM-Newton needs only three absorbers to be fitted. 
Changes of ionization state in at least one absorbing zone are re- 
quired to fit data taken at different epochs. The partial-covering ab- 
sorption model might offer an alternative explanation for the vari- 
ability and spectrum of MCG-6-30-15, but the reflection model, 
together with the light-bending model, is fully consistent with the 
data and provides a more physical point of view. 



2 DATA REDUCTION 

Data from XMM-Newton, Suzaku and Chandra provide observa- 
tions of different epochs of MCG-6-30-15. These observations are 
in different flux states so that we can study spectrum variations and 
evolution in recent years. The 2-10 keV bolometric flux is 4.26 x 
10"" ergs cm"^ s"^ for the XMM-Newton PN, 4.14 x 10"" ergs 
cm~^ s"^ for the Suzaku FI XIS, and 4.03 x 10"" ergs cm"^ s"^ 
for Chandra HEG. The simultaneous BeppoSAX observation with 
XMM-Newton and PIN data of Suzaku offer hard X-ray data above 
10 keV and we use them for the study of hard excess. 



2.1 XMM-Newton and BeppoSAX 

We use the XMM-Newton observation taken from 2001 July 31 to 
August 05. In this observation EPIC MOS and PN cameras were 
operated in small window mode. The observation data files (ODFs) 
were reprocessed with the standard procedures of XMM-Newton 
Science Analysis System (SAS V9.0.0). Source data were extracted 
by XSELECT using a circle with a radius of 30 arcsec. The HEASoft 
package is incorporated into SAS. When saving the extracted spec- 
trum in XSELECT, the script xsl_xmin_epic_makeresp runs 
and creates the redistribution matrix file (rmf) and the ancillary re- 
sponse file (are). Background was taken from a region within the 
small window using a circle of the same size. Using the epatplot 
tool in SAS we confirm the observation is not piled up. Total expo- 
sure time of the PN camera is about 227 ks. The effective energy 
band of EPIC cameras is 0.5-10 keV. 

We reduced the Reflection Grating Spectrometer (RGS; 0.4- 
2.0 keV) data following the standard procedures. Both the first or- 
der spectra of RGS 1 and RGS2 data were used. Data of these in- 
struments overlap over most of their effective energies. To analyse 
the data more efficiently we used RGS 1 spectrum between 0.4 and 
1.7 keV, and RGS2 spectrum between 0.9 and 1.2 keV. We use only 
RGS data, which are with better energy resolution, at low energies 
in all our fittings. 

To study the spectral variability, difference spectra were pro- 
duced by subtracting the low-flux state spectrum from the high- 
flux_state one. Flux cuj.s_were chosen in a similar way mentioned 
in fVaughan & FabianI ( 12004 ). We divided the EPIC PN 0.2-10 keV 
light curve into two flux intervals by 20 counts s^^ so that the two 
slices have a comparable number of total counts (~ 3 x lO''). The 
slice with the count rate lower than 20 was used as the low-flux 
state spectrum, while the upper slice is devoted the high-flux state 
spectrum. 

BeppoSAX observed the source during the 2001 XMM- 
Newton observati on. Data were reduced following the guidelines in 
iMatt et all ( Il997h . We use data from the Phos witch Detector Sys- 
tem (PDS; 13-200 keV) to extend the coverage of the observation 
to the higher energy band. The total exposure time of the PDS data 
is about 49 ks. 



2.2 Suzaku 

In 2006 January, MCG-6-30-15 was observed three times by 
Suzaku for a total of about 338 ks of exposure time. The X-ray 
Imaging Spectrometer (XIS) was operated in the normal mode. 
5x5 and 3x3 editing modes were used throughout the observation. 
We use HEASoft V6.8 software package provided by NASA to reduce 
the data following the Suzaku Data Reduction Guide. We created 
new cleaned event files using the XISPI tool. Source products 
were extracted by xselect using a circular region of 4'. 3 radius, 
while background products were extracted from a smaller circular 
region away from the source. Response files were produced by the 
script XISRESP, which calls XISRMFGEN and XISSIMARFGEN 
automatically. Spectra of the three front-illuminated (FI) CCD XIS 
detectors (XISO, XIS2 and XIS3) were combined using the ftool 
addascaspec. For the difference spectrum, the low-flux state 
and high-flux state were defined in the same way as for XMM (see 
Fig.E). 

The Hard X-ray Detector (HXD) was operated in XIS- 
nominal pointing mode. The background spectrum of HXD/PIN 
consists of a non-X-ray background and a cosmic X-ray back- 
ground. We used the non-X-ray background event file from a 
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Figure 1. The light curves on this figure ai'e of the three satellites we use. 
The top panel is the XMM-Newton PN 0.2-10 keV light curve, while the 
middle one belongs to the Suzaku FI XIS detectors (XISO, XIS2 and XIS3) 
of 0.5-12 keV, and the bottom panel shows the combined Chandra 0.45-10 
keV light curve of first order HEG and MEG instruments. The dash line 
in each panel hes on the edge of low-flux and high-flux states. The value 
was chosen to make the low-flux and high-flux slices have comparable total 
number of photons (~ 3 X 10^ fox XMM-Newton, ~ 1.6 X 10^ for Suzaku, 
and ~ 2.6 X lO'^ for Chandra). 



database of background observations made by the PIN diode. The 
appropriate response files for the observation epoch were used as 
well. The cosmic X-ray background was obtained by model esti- 
mation using the PIN response for flat emission distribution. We 
combined the source spectra and non-X-ray background spectra 
of different datasets first, then created a cosmic X-ray background 
spectrum for the entire observation, and added it to the total non- 
X-ray background spectrum. 

2.3 Chandra 

The Chandra HETGS observed MCG-6-30-15 between 2004 
May 19 and 27. Tools in the latest CIAO V4.2 software package 
were used to reduce the data. We started with level 1 raw data 
and followed standard procedures in the analysis guide Q. We 
examined the background light curve and removed periods affected 
by background flaring, and obtained a good exposure time of about 
497 ks in total. The tgdetect tool was used to find the zeroth 
and ±1 order position. With appropriate filters, a level 2 event 
file could be generated and a type ll pha (hereafter pha2) file 
which contains complete information could be extracted. Different 
orders of spectra of both the medium-energy gratings (MEG) 
and high-energy gratings (HEG) were separated from the pha2 
file. In order to do the grating analysis in XSPEC, we also split 
the background file from the pha2 file by executing the tg_bkg 
script. RMF and ARF files were created by the mkgrmf and the 
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Figure 2. The figure shows the difference spectra of the three satellites 
plotted as energy against the data/model ratio. The model of a power-law 
with Galactic absoiption modification fitted across 3-10 keV (3-5 keV for 
Chandra MEG). The difference spectra from top to bottom belong to XMM- 
Newton PN camera, Suzaku FI XIS detectors and Chandra HETGS, respec- 
tively. The absoiption structure around 1.8 keV in the middle panel is pos- 
sibly caused by the detector These difference spectra show a highly siirtilar 
structure between 0.5 keV and 5 keV. 



fullgarf tools, respectively. Since the zeroth-order source is 
piled up, only ±1 order spectra were used for further study. The 
combined HEG (0.8-7.5 keV) and MEG (0.4-5 keV) first-order 
light curve is shown in Fig.[T] 



3 SPECTRAL FITTING 

3.1 A Quick Look at Difference Spectra 

We create a new spectrum by subtracting a low-flux spectrum from 
a high-flux spectrum, where the definition of low-flux and high- 
flux states has been elaborated in Section|2] The resulting spectrum, 
which is called a difference spectrum, offers information on the na- 
ture of the variability because constant components in the original 
spectrum are removed. The difference spectrum was initially intro- 
duced into the analysis of MCG-6-30-15 as a demonstration that 
the variability is indeed consisted with a pure powerlaw (as further 
shown here). The absoiption is a multiplicative component so in the 
case of a variable powerlaw, absorption components are revealed 
clearly in the difference spectrum (Fig.|2)- If the variability is more 
complex, and no simple spectrum emerges , then a diff erence spec- 
trum is difficult to int erpret. Both Fabian et al.l (l2002l) (using data 
of XMM-Newton) and lMiniutti et all ( l2007h (using data of Suzakit) 
reported that the reflection component remains constant through- 
out the observation period and does not appear in the difference 
spectra. In addition, the energy band above 3 keV in the difference 
spectra can be simply fitted by a power-law with Galactic absorp- 
tion (Nn = 4.06 X 10^"cm"^), which means the warm absorbers 
have little effect on the spectrum above 3 keV, particularly in the 
Fe-K band, in contrast with the band below 3 keV. 
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Figure 3. The black data points in the upper panel belong to the MEG 
spectrum, while the red ones are of the HEG spectrum, and the blue line 
stands for the best-fitting model, which consists of a powerlaw continuum, 
a ionised and a neutral reflection components, an additional naiTow Gaus- 
sian line, and four absorbing zones. The lower panel shows the data/model 
ratio of this fitting, and the data points have been re-binned for clarity. De- 
tailed information of this fitting is listed in Table [T] and |2] The model is in 
excellent agreement with the data above 1 keV and struggles slightly below 
1 keV, but main features are fitted. The small box shows the zoom-in iron 
line band. 

We reproduce the difference spectra of XMM-Newton and 
Suzakii and confirm the conclusions in previous work. The differ- 
ence spectra of Chandra HETGS are generated to compare changes 
in the difference spectra from observations taken at difference 
epochs. The difference spectra of the tliree satellites are shown as 
a ratio to a simple power-law model with Galactic absorption in 
Fig-IH They are remarkably similar at low energies; they are much 
alike below 3 keV, at least down to 0.5 keV, which is the lowest 
energy we can trust XMM-Newton and Suzaku. An obvious drop 
below 2 keV is seen and goes to nearly the same depth in all differ- 
ence spectra. This hints that the low energy spectra we obtained can 
be interpreted by the same model, and the warm absorbers did not 
change much between epochs. The 3-10 keV (3-7.5 keV for Chan- 
dra) difference spectra can be fitted by a simple powerlaw, imply- 
ing the warm absorbers cause little curvature above 3 keV. There is 
no clue that warm absorbers mimic the "red wing" structure at iron 
line energies. 

3.2 Warm Absorbers 

The spectrum of MCG- 6-30-15 s h ows c omplicated absorption 
in the soft X-ray band. iLee et aL I ( l200l[) showed at least two 
ionisation zones are needed to interpret the absorption lines. 
ISako et al] ( l2003h verified two outflow components with veloc- 
ities of -150 km s~^ and -1900 km s~^. Later studies (e.g. 
iTumer et al ]|2004lYoung etal.l2005ll also confirmed the existence 
of the two kinematic components, with slightly different veloci- 
ties. McKeman et al. (2007) reported both com ponents and found 
two ionisation states in the slow components. IMiller et alJ ( 1200 Sh 
adopted the three absorption zones mentioned in previous literature 
and added two more offset comp onents to yield a goo d fit. The lat- 
est detailed kinematic analysis of Holczer et al. ( 2010), which used 
the absorption measure distribution (AMD) method which enables 
a full reconstruction, indicated there is a "local" {z = 0) component 
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Figure 4. The upper figure shows the best fit of the simultaneous XMM 
and BeppoSAX observation, and the lower one of the Suzaku observation. 
The black data points of XMM in the upper figure are heavily re-binned for 
clarity, and data points of BeppoSAX are shown in red. However, the red 
data points in the lower figures are of Suzaku XISl detector, whose CCD 
is back-illuminated and collects better data at lower energies than those 
of front-illuminated detectors. The high energy spectrum (above 10 keV) 
of Suzaku PIN is shown in green in the lower figure. The blue line again 
represents the theoretical model. It could be clearly seen that the model fits 
the high energy data well in both figures. The small boxes inside the figures 
as well show the zoom-in iron line profile. 



with a column density comparable to that of the Galactic absorp- 
tion column of this object. The newly found third kinetic system 
(~-2300 km s-\ the same shift as the redshift z = 0.007749 of 
MCG-6-30-15) was confused with the fast component (~ -1900 
km s ~^) because o f similar velocities in the local frame of refer- 
ence. iHolczer et al. 

I UolO) explained the component as absorption 
by ionised Interstellar Medium (ISM) in the Milky Way or areas 
nearby. So far the local absorption in MCG-6-30-15 is not yet con- 
firmed by other work and its turbulence velocity remains unknown. 
Summarising the previous work, a highly-ionised (log ^ > 3.5) 
absorption zone with turbulence velocity utmb = 500 km and 
at least two zones of Wtu* = 100 km s~^ with different ionisation 
states are required to explain t he Chandra HETGS data (see Table 
5 and 6 in lHolczer et alj20I(ll) . 

We generated the ionised absorber models using xstar V2.2.0. 
The source luminosi ty between I and 1000 Ryd was assumed to 
be 2 X 10*^ erg s~^ dYoung et alj2005b with a powerlaw spectrum 
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Figure 5. These figures show the zoom-in iron line fitted by the best-fitting model. Line profiles from left to right belong to Chandra HEG , XMM PN, and 
Suzaku respectively. 



shape at 0.1-20 keV, the gas density was set to be 10^^ cm~'^, and 
the covering factor equal to 1, for all models. Fitting parameters 
in these models are the photon index, the absorbing column, the 
ionisation parameter log ^, and the iron and oxygen abundances 
except for the fast component, whose oxygen abundance was fixed 
to unity. Previous work adopted a photon index F between 2.0 and 
2.2 when calculating the absorption tables. We allow F of absorbers 
to vary with the value of the powerlaw incident continuum. Iron and 
oxygen abundances are not fixed in case of super-solar or sub- solar 
abundant environment. Turbulent velocities titmt = 100, 500, and 
1000 km s"^ and temperatures T = 10"*, 3 x 10*, and 10^ K were 
applied to produce different absorption models that cover a wide 
range of parameters. 

We then tested these absorbers on the Chandra MEG spec- 
trum, with three absorbing zones and a blurred reflection compo- 
nent as a start. Since the outflow velocities of both the slow and the 
fast comp onents determined b y previous work are in good agree- 
ment (see lHolczer et alj|2010h . we constrained the redshifts of the 
absorbing zones to vary in a small range. The primary fitting results 
indicate a zone with Wmrb = 500 km s^^ and two with UtuA = 100 km 
s~^ are required to fit the spectrum, which is consistent with previ- 
ous work. The model was further tested on the low energy spectra 
(< 3 keV) of XMM and Suzaku, and we found the ionisation state of 
the fast component is nearly constant on all dat asets. Therefore we 
fixed the ionisation parameter log ^ to be 3.82 jMiller et alj|2007l : 
iHolczer et al ]l2O10l) in all the following fits. 

A fourth absorbing zone, corresponding to the third kinetic 
system detected in Holczer et al. (2010), with utuA = 100 km s^^, 
2 = 0, and A^'h = 4.06 x lO'^" cm~^, was added to test how much 
this component affects the fitting, and the detailed information of 
the absorbing zones we use is listed in Table[T] The results of other 
fitting parameters are shown in Table |2] Clearly the model with the 
fourth component yields a better fit to the Chandra spectrum of 
MCG-6-30-15. The model struggles slightly below 1 keV, where 
the most residuals come from. We then model the Fe L edge at 0.7 
keV using dust in tbnew, and the result is shown in Fig. [3] 

As to the difference spectra, we fit them with a model consist- 
ing of a powerlaw. Galactic absorption and the absorbing zones we 
use to fit the real spectrum. The model fits all our difference spec- 
tra without difficulty, indicating the reflection component is suffi- 
ciently removed by the difference spectrum technique. 



3.3 Full Spectrum 

When constructing a model for the full spectra, a relativistic blurred 
reflection component and a neutral reflection component should be 
included in order to fit the skewed, broad iron li ne and a narrow iro n 
line reflected from distant matter (detected by lYoung et al.ll2005l) . 
We model reflection from the innermost re gions of the accretion 
disc with the self-consistent model reflionx jRoss & Fabianl2005h 
convolved by the kdblur model, and the neutral reflection with re- 
flionx only = 1.0). The iron abundance Ape of the reflection 
components are bound for consistency. We set the outer radius of 
the accretion disc to be 400 Rg in all fits, where Rg is the gravita- 
tional radius GM/c^. Nevertheless, the two-reflection model is still 
not enough to fit the iron line structure perfectly, and we introduce 
another narrow line modelled by ZGAUSS. The new model works 
well on joint MEG and HEG spectra (see Fig. [3} and explains the 
complex structure of the iron line. 

Some claimed the hard excess seen in some AGNs, included 
MCG-6-30-15, cannot be explained by th e reflection model , and 
only clumpy absorbers model can work. Miller et all (|2008|) first 
modelled the hard exc ess with absorbed reflection, but in their later 
work (i.e. lMiller et al.li2 009) they claimed the high energy emission 
is dominated by absorption, and the hard excess is reproduced by 
the continuum covered by a partially covering absorber To check 
the availability of the reflection model at high energies, we use the 
simultaneous XMM and BeppoSAX spectra and both Suzaku XIS 
and PIN data, adopting the same model used to fit the Chandra 
HETGS spectra. The results (see Fig.|4) show the reflection model 
is able to produce a "Compton hump" easily, and the hard excess 
can be well explained by simple reflection without any additional 
absorption. The zoom-in iron lines of different datasets are shown 
in Fig.|5] and Fig.|6]shows the decomposed components of the final 
model. 

The resulting xt shown in Table |2] may not be satisfactory 
enough. If 1 percent systematic errors are included (using the sys- 
tematic command in xspec) for the whole band of Suzaku, and for 
XMM and Chandra above 1 keV, xl is acceptable at 1.00, 1.02 and 
1 .04, respectively. 



3.4 More Tests on XMM and Suzaku 

We test the 3-absorber model on XMM and Suzaku observations. 
Since data points at low energies are statistically stronger than those 
at higher energies, the model tends to fit the low-energy spectrum 
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Table 1. The table lists parameters of the warm absorbers in the best-fitting 
model of different satellites, in which A^h is given in 10^^ cm^"^, temper- 
ature in Kelvin, and ? in erg cm s^^. The column density of the local ab- 
sorption co mponent is fixed at th e same value with Galactic absorption, as 
suggested in lHolczer et al.l ilOlOl) . in all fittings. Meanwhile, the ionisation 
state of the fast component is fixed as well. 



Absorber 


fast 


slow (1) 


slow (2) 
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27 
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1 61+°°^ 


"•^^-0.14 
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and results in an unfitted iron line profile in Suzaku data. We then 
fit the spectrum above 3 keV first and fix the parameters of the con- 
volution model (the emissivity profile index, inner radius Rm and 
inclination angle <f)) so that the model will not ignore the iron line 
completely when fitting the full spectrum. The results are shown in 
Tabled 

As stated in Section [T2l Chandra HETGS data tend to include 
the local absorption component, and we test how necessary this 
component is needed on XMM and Suzaku observations as well. 
The of the fittings of both XMM Suzaku data do improve (see 
Table |2]l, though not much in XMM data. However, the ionisation 
state of one of the slow sub-components in Suzaku is close to the 
value of the local component (see Table [T), so it remains unclear 
whether Suzaku data are able to resolve the additional absorption. 
Whether the 4-absorber model or the 3-absorber model is used, the 
ionisation states of the slow component remain constant, indicating 
the local absorption may not much affect the fits. 

The joint XMM and BeppoSAX and the Suzaku observations 
do not show a strong tendency for the local absorption component. 
To test if it is likely th at even fewer absorbing zones, say, two as 
iMivakawa et^ j2009 f) suggested, explain the spectrum, one of the 
slow sub-components was removed from the model. The results of 
the 2-absorption model are also shown in Table|2] The 2-absorption 
model results in a worse in both XMM and Suszaku datasets, im- 
plying at least three absorbing zones are required to fit the spectrum 
of MCG-6-30-15. 

To examine if the real spectra in different periods have simi- 
larities like the difference spectra, we put the best-fitting model of 
the joint XMM and BeppoSAX spectrum to the Suzaku broadband 
spectrum and modify the normalisation slightly. It seems the spec- 
tra of these observations above 3 keV look fairly similar, but there 
are obvious differences below 2 keV, where the absorption domi- 
nates the spectrum. The fact that the iron line band and hard excess 
appear constant but variations in source luminosity b elow 10 keV 
motivated the work of the clumpy absorption model. iMiller et"ail 
( 12009 ') suggested the long-term X-ray variability may be caused by 
variation in the covering factor. 

Nevertheless, the differences below 2 keV in the spectra may 



Table 2. The Table below lists the fitting parameters and obtained by 
the reflection model with different numbers of absorbing zones. The top 
block shows the result using a reflection model with four absorbers; the 
middle block displays the result fitted with a 3-absorber model (local ab- 
sorption taken out from the 4-absorber model); the bottom block shows the 
result using a reflection model with one slow component and one fast com- 
ponent only. In the case of 2 warm absorbers, only XMM and Suzaku data 
have been tested. As mentioned in Section [T4l emissivity profile index, iron 
abundance and inclination angle are fixed when testing XMM and Suzaku 
observations. 



parameter 


Chandra HETGS 


XMM + BeppoSAX 


Suzaku 


model with 4 absorbing zones 


r 

index 

fii„(Rg) 


1.97 ±0.00 
8.00+n^ 

1 qi+0.08 

I'-^l-O.OO 
1 70+0-10 

fi ^is+o-oe 

■J3->J_35.0 


^■'^'-'-0.01 
q 70+0.05 
•^•'"-0.08 
-, C7+0.13 
l-°'-1.57 
-1 7q+0.19 
'•^-0.12 

6.52 ±0.03 
37.7l^;*° 


1.98 ±0.01 

3.09 

2.50 
4 00+"-°° 

f. qO + 0.01 

44.0° 


xVd.o.f. 


2417.7/2139 


5059.3/3809 


1684.7/1576 


model with 3 absorbing zones 



r 


1 94+0-01 
^■^^-0.00 


1.96 ±0.01 


1.98 ±0.01 


index 


s.ooins 


3 84+0-09 


3.34 


fiin(Rg) 


T qT+0.04 
^■■'^-0.00 


1 O1+0.41 
1-'^1-1.31 


2.66 


Ape 


2.07«:1| 


2 77+0.79 
^- "-0.65 


4 00+"-°° 
^•1^-0.11 




6.48t°:°* 


"■^^-0.03 


f. qO + 0.01 





qc 7+1.4 




41.6° 


xVd.o.f. 


2516.4/2140 


5111.9/3810 


1688.3/1577 



model with 2 absorbing zones 



r 

index 
Rin(Rg) 

xVd.o.f. 



1.98 ±0.01 
4 40+°-i2 

^+0.33 
-1.24 

r| + 0.01 

-0.00 
°-^*-0.05 

6040.9/3815 



1.24: 
1.00"* 



1.98 ±0.00 
2.98 
3.11 
4 00+"-°° 

^•'-"J-0.17 

6.37 ±0.02 

33.0° 
1692.0/1582 



be alternatively ex plained by changes in column densities (see 
I Walton et al1l2010h and ionisation states in warm absorbers. Since 
there are outflowing components in MCG-6-30-15, column densi- 
ties of these components may vary throughout time, and their ioni- 
sation states may possibly change due to the re-allocation of clouds. 
We fit Suzaku data by the best-fitting model obtained from the joint 
XMM and BeppoSAX spectrum, with all the parameters of warm 
absorbers fixed except the column densities. 



4 DISCUSSION 
4.1 Warm Absorbers 

Table |2] shows the of the fits using the model with and without 
the fourth local absorption component. Only Chandra data show a 
strong tendency to include the local absorption component, which 
improves the fitting of Suzaku data but does not show strong evi- 
dence of existence, and XMM observation only gets a slightly bet- 
ter fit with it. The local absorption component in fact contributes 
little to the spectrum, and only a few oxygen lines around 20 A 
are slightly shifted by the component. Besides, if this component 
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Table 3. The tabl e shows a comparison among different models. 
iHolczer et al and present work use full covering zones only, but 

the former used a phenomenological cubic spline conti nuum, and the later 
models the continuum with a powerlaw and a reflection. iMiller et alj j2009[) 
used both full and partial covering absorbers, and a simple continuum with- 
out reflection. A^h is given in 10^^ cm^'^, temperature in Kelvin, and ? in 
erg cm s^^. 



work 




Holczer et al. 


Miller et al. 


present work 


full covering zones 


fast 
slow(l) 
slow(2) 

local 


Nn 
log? 
Nh 
log? 

Nh 

log? 

Nh 

log? 


81 ± 7 
3.82 ±0.03 
2.3 ±0.3 
-1.5-0.5 
3.0 ±0.4 
1.5-3.5 
0.40 


(80.0) 

(3.95) 
0.27 ±0.03 
0.88 ±0.16 
11.8 ±0.5 
2.39 ±0.01 


209 4+36.9 

(3.82) 
40+0.31 
'^•^'-'-0.42 

1.71 ±0.03 

n. 97+0.20 
U.^(_0.i3 
2 47+0.03 
^■^'-0.16 

(0.406) 

-1 62+°-°"'' 


partial covering zones 


zone 1 
zone 2 


log? 

Nh 
log? 




1910 ± 300 

29 ±1 
1.38 ±0.03 





is caused by the ISM in the local group as iHoIczer et all ilOUj) 
suggest, its ionisation state should stay constant in all datasets. But 
the ionisation states of the component obtained from Chandra and 
Suzaku are significantly different (see Table [T), indicating the data 
may have trouble constraining the parameters of this zone. 

The result is expected because XMM and Suzaku do no t 
have the sharp energy resolution of Chandra. lHolczeretai] ( l20ld) 
claimed the third kinetic component was previously confused with 
the fast component because of similar velocities; even Chandra 
does not have enough resolution to constrain the line width. There- 
fore it is not surprising that Suzaku could not discern these two 
components, considering the weak signature one of them. Our spec- 
tral fitting results do favour its existence in the Chandra spectra of 
MCG-6-30-15, but cannot judge whether the component is caused 
by the Galactic absorption or another outflow; the issue is not rel- 
evant to the present work. So far there have already been three ki- 
netic systems, which are described by four ionisation absorbers in 
this work, confirmed to be responsible for the soft X-ray band. The 
highly ionised fast component, with umrb = 500 km s^^, is re- 
quired in all the observations, though the column density needed 
for each observation differs (see Table [T). The fast component was 
reported in all previous work of kinetic analysis of MCG-6-30-15, 
and our result again shows its importance in modelling the spectra. 
The main features it causes are the Fe xxv (~ 6.7 keV) and Fe xxvi 
(~6.9 keV) absorption lines, which can be seen in Fig.|5] The fast 
component is also the most important absorber that affects the iron 
line band. 

Summarising the fitting results obtained by the 3-absorption 
and 2-absorption models (see Table|2}, it seems XMM and Suzaku 
data are able to distinguish the sub-components in the slow outflow. 
Nonetheless, the values of ionisation states in XMM and Suzaku 
data are different from that in Chandra data (see Table [T]!. This 
could be because changes in ionisation states are needed to inter- 
pret the variation in the low spectrum, or simply because XMM and 



Mode 




1 10 100 

Energy(keV) 



Figure 6. The figure shows the decomposed components of the 4-absorption 
model. The black dot-dash line above shows the reflection component con- 
sisted of a blurred reflection, a neutral reflection and a narrow Gaussian Une. 
The blue dash line represents a F = 2.0 powerlaw and the red solid line 
shows the effects of the warm absorbers. Note that the reflection component 
has been moved above for clarity. 



Suzaku spectra are not good enough to constrain parameters of the 
sub-components well, though they do need two ionised zones to 
model the slow component. 



4.2 Robustness of the Reflection Model 

The p revious work (e.g. IVaughan & Fabiaiil 120041 ; iMiniutti et al.l 
I2OO7I etc.) which applied a reflection model on MCG-6-30-15 sug- 
gested the photon index to be F 2 — 2.2; the emissivity index is 
around 3; the inner radius of the disc lies at somewhere less than 
2 Rg; an iron abundance of 2 to 3 solar is needed to interpret the 
spectra, and the inclination angle should be around 40°. Most of 
the fitting parameters of the relativistic reflection model, shown in 
Table [2I are in good agreement with the assumption of previous 
work. Note that the previous work analysed the spectra above 3 keV 
only, and our work includes the low energy band that is seriously 
influenced by the warm absorbers. This implies the results obtained 
from the reflection model would not be affected by absorbing zones 

and are reliable. 

IZvcki et alJ 1 I2OIOI) claimed the results received through the 
reflection model are model-dependent. They found quite different 
powerlaw indexes when they fit the Suzaku spectrum above 3 keV 
(F > 2.1) and the spectrum above 1 keV (F < 2.0), but we 
do not find this with our analysis. The F we obtained from a 3- 
absorber model or a 4-absorber model are all between 1.9 and 2.1, 
and the energy band we use is down to 0.5 keV. However, photon 
index F does change slightly when we adopt models with fewer 
warm absorbers. The disagreement might be caused by the warm 
absorbers u sed to model the so ft X-ray spectrum. The diverse pho- 
ton indexes IZvcki et al.l ilOVj) obtain is simply due to the warm 
absorbers rather than the reflionx model. Additionally, the extra 
neutral (^ = 1.0) reflionx component that we used to fit the nar- 
row iron line works well, implying reflionx can produce neutral 
or mild ionised reflection without difficulty (we use a reflionx 
grid w hich extends to ^ = 1.0, contrary to that used bv lZvcki et al.l 
I20T0I) . 
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4.3 Comparisons with Previous Research 

Table [Sllists parameters of warm absorbers obtained by different 
models. [Holczer et al. (2010) suggested the slow component spans 
a considerable range of ionisation, which are -1.5 < log ^ < 0.5 and 
1.5 < log ^ < 3.5 for the sub-components, re spectively. The ioni- 
sation state of slow (1) in lMiller et al.l ( |2009|) and our best-fitting 
model for Chandra data (see Table O do not fall in the former 
range. The ionisation states of slow (2) of Miller et al. (2009) and 
present work are qu ite close and lie well in the expected range of 
lHolczeretal] ( l2010l) . 

The column density we get for each absorbing zone differs sig- 
nificantly from the values obtained in previous literature. They are 
probably the most difficult parameters to be constrained, because 
results depend on the way of modelling the continuum and the re- 
search tools that have been used. Column densities are also some- 
times under- or over-estimated in order to fit complex structures 
in a spectrum as com plicated as in MCG-6-30-15, as discussed by 
iHolczer et al.l ( |2010|) . The different parameters we used for gener- 
ating XSTAR grids might again be the reason we have disagreement 
with previous estimations. Moreover, reflionx contributes some 
soft emission and warm absorbers may need higher column densi- 
ties to fit the spectrum. 

iMiller etai] ( I2008ll2009h used two additional partial-covering 
clumpy absorbing zones to model the spectrum. One of them mim- 
icked the shape of broadened asymmetric iron line, and the other 
partially covered the continuum in order to explain the hard excess. 
The main difference between the absoiption-dominated model and 
the reflection-dominated model is that t h e form er has no distortions 
due to relativistic effects. lMilleretal.l ( l2009l) claimed the reflec- 
tion (relativistic effects not considered) model failed to interpret 
the hard excess. In fact the hard excess can be simply explained by 
a relativistically blurred reflection component without any absorp- 
tion. The complex iron line structure can also be fitted by a combi- 
nation of blurred and neutral reflections and absorptions caused by 
a highly ionised warm absorber. Both the absorption and the reflec- 
tion models explain the spectrum and variability of MCG-6-30-15. 
Our result does not rule out the "3-1-2" model but shows the reflec- 
tion model that includes relativistic effects works. We consider this 
to be a more consistent physical scenario. Since the energy release 
in both models occurs very close to the black hole, where relativis- 
tic reflection is expected. 



5 CONCLUSION 

The complex X-ray spectrum of MCG-6-30-15 can be well ex- 
plained by a model consisting of relativistic reflection and three or 
four comp letely covered abso r bing z ones, in contrast with the con- 
clusion of iMiller et al. I ( l2008ll2009l) . Our result does not rule out 
the absorption-dominated scenario but does support the reflection- 
dominated model. In addition, we use absorbing zones alread y con- 
firmed by previous analysis, while IMiller et al.l(l2008Ll2009h added 
two offset components to yield good fits. The "3-h2" model is sup - 
ported by the Principal Components Analysis dMiller et al.ll2007l) . 
However, if there is an absorber responsible for the variability and 
the "red wing" structure at 2-6 keV, structures such as curvature 
should appear in the difference spectra at this band, and it should 
leave some signature that could be traced by kinetic analysis. The 
difference spectra retrieved from the observations we use (see Fig. 
|2]l is well fitted by a simple powerlaw above 3 keV, contrary to the 
assumption that the spectrum above 3 keV is dominated by absorp- 
tion. 



A highly-ionised fast outflow is confirmed by all observations. 
This is the most constrained absorbing zone among the warm ab- 
sorbers needed to explain the X-ray spectrum of MCG-6-30-15. 
The slow component with two ionisation states are confirmed as 
well, but XMM and Suzaku data are not good enough to determine 
the ionisati on states correctly. Th e "local absorption" component 
detected by iHolczer et all ( 1201 Oi) does improve our fitting to the 
Chandra HETGS spectra. By summarising the results we find three 
absorbing zones are required to fit XMM and Suzaku data, and a 
fourth zone is need t o fit Chandra data com pletely. The 2-absorber 
model suggested bv lMivakawa et al" I (2009^ is not sufficient to in- 
terpret the complex spectrum of MCG-6-30-15, so the discrepancy 
in photon index between the spec trum above 3 keV and above 1 
keV claimed bv lZvckietal] ( l201(]h is due to the use of inappropri- 
ate warm absorbers. 

We also test the robustness of our model and find the results 
are in good agreement with previous work. Therefore we conclude 
that the reflection model can robustly interpret the spectra of MCG- 
6-30- 1 5 in the range from 0.5-200 keV without any partial-covering 
absorbers. 
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